Measurement of quasi-elastic 12 C(p,2p) scattering at high momentum 

transfer 

Y. Mardor 1 , J. Aclander 1 , J. Alster 1 , D. Barton 2 , G. Bunce 2 , A. Carroll 2 , N. Christensen 3 *, H. Courant 3 , S. 
Durrant 2 , S. Gushue 2 , S. Heppelmann 4 , E. Kosonovsky 1 , I. Mardor 1 , M. Marshak 3 , Y. Makdisi 2 , E.D. Minor 4t , I. 
Navon 1 , H. Nicholson 5 , E. Piasetzky 1 , T. Roser 2 , J. Russell 6 , C.S. Sutton 5 , M. Tanaka 2 *, C. White 3 , J-Y Wu 4§ . 



O 



> 

O 
O 

o 

o> 
">< 

=5 



X 



1 School of Physics and Astronomy, Sackler Faculty of 
Exact Sciences, Tel Aviv University. 2 Brookhaven Na- 
tional Laboratory. 3 Physics Department, University of 
Minnesota. 4 Physics Department, Pennsylvania State 
University. 5 Mount Holyoke College. 6 Physics Depart- 
ment, University of Massachusetts Dartmouth. 

We measured the high-momentum quasi-elastic 12 C(p,2p) 
reaction (at 6» cm ~ 90° CM.) for 6 and 7.5 GeV/c incident 
protons. The three-momentum components of both final state 
protons were measured and the missing energy and momen- 
tum of the target proton in the nucleus were determined. 

The validity of the quasi-elastic picture was verified up to 
Fermi momenta of about 450 MeV/c where it might be ques- 
tionable. Transverse and longitudinal Fermi momentum dis- 
tributions of the target proton were measured and compared 
to independent particle models which do not reproduce the 
large momentum tails. We also observed that the transverse 
Fermi distribution gets wider as the longitudinal component 
increases in the beam direction, in contrast to a simple Fermi 
gas model. 



Quasi-elastic (QE) scattering is a process in which a 
projectile is elastically scattered from a single bound nu- 
cleon in the nucleus, which we call the "target nucleoli", 
while the rest of the nucleus acts as a spectator. Specifi- 
cally, the QE (p,2p) scattering at large momentum trans- 
fer provides a method for measuring the high momentum 
tails of the nuclear wave function. This fact can be under- 
stood by considering the s-scaling law for high momen- 
tum transfer hadronic reactions M. The elementary pp 
clastic differential cross section scales as da/dt ~ 1/s 10 
for fixed (s/t), where s and t are the Mandelstam vari- 
ables. Farrar et al. Q pointed out that this scaling will 
cause the QE (p,2p) reaction in the nucleus to favor 
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strongly the scattering from those target nucleons that 
are moving in the beam direction with large Fermi mo- 
mentum , because the value of s is reduced under those 
kinematic conditions. The same reaction also attracted 
much attention in recent years Q| in connection with the 
QCD prediction that the nuclear attenuation will van- 
ish at asymptotically large momentum transfer, a phe- 
nomenon called "Color Transparency" 0]. Thus, at large 
but finite energies the QE (p,2p) reaction , with emerging 
protons of a few GeV/c, is a good tool for experimental 
studies of the nuclear transparency. Those results will be 
discussed in a forthcoming publication || . 

This letter will concentrate on the identification of QE 
events and will test the validity of the QE picture up 
to large Fermi momenta, where it might become ques- 
tionable. If the QE picture is valid, then over the mea- 
sured kinematical range it is possible to separate nuclear 
properties from the nuclear reaction mechanism. We will 
also present transverse and longitudinal Fermi momen- 
tum distributions of the target proton and compare them 
to independent particle models. 

We measured the high-momentum transfer quasi- 
elastic (p,2p) reaction at 9 CTn ~ 90° on carbon for 6 and 
7.5 GeV/c incident protons in a kinematically complete 
coincidence experiment. The three- momentum compo- 
nents of both high p t final state protons were measured, 
which yielded the missing energy and momentum of the 
target proton in the nucleus. 

The experiment (E850) was performed at the AGS ac- 
celerator at Brookhaven National Laboratory with the 
EVA spectrometer The spectrometer consists of 

a super-conducting solenoidal magnet operated at 0.8 
Tesla. The beam enters along the z axis and hits a series 
of targets located at various z positions. The scattered 
particles are tracked by four cylindrical chambers (Cl- 
C4). The radii of the cylinders range from 10 to 180 
cm. All cylinders and targets can be moved along the 
solenoid axis in order to optimize the angular acceptance 
range for each beam momentum. Cylinders C2-C4 have 
4 layers of 2 m long straw drift tubes, whose diameters 
range from 1 cm for C2 to 2 cm for C4. The high resis- 
tance central wires are read out at both ends, providing 
position information along the z direction. Thus, one can 
extract the z position of the particles in the cylinders as 
well as their azimuthal angles as they are bent in the 
axial magnetic field. This provides the transverse mo- 
mentum of the particles and their scattering angle. The 
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1 m long CI cylinder with a tube diameter of 0.5 cm, was 
read out at one end only. The readout electronics were 
designed especially for this spectrometer |J. The straw 
tubes were filled with a 50/50 mixture of argon-ethane 
gas at atmospheric pressure. The drift time measure- 
ment from the central wire had a spatial resolution of 
about 0.3 mm. Three solid targets, CH 2 , C and CD 2 
(enriched to 95%) were placed on the z axis inside the 
CI cylinder separated by about 20 cm. They were 5.1x5.1 
cm 2 squares and 6.6cm long in the z direction except for 
the CD 2 target which was 4.9 cm long. Their positions 
were interchanged at several intervals in order to reduce 
systematic uncertainties and to maximize the acceptance 
range for each target. 

The spectrometer was located on the secondary line 
CI of the AGS. The beam passed through a sequence of 
two differential Cerenkov counters which identified the 
incident particles. The beams ranged in intensity from 1 
to 2 • 10 7 over a one second spill every 3 seconds. Two 
counter hodoscopes in the beam provided beam align- 
ment and a timing reference. Three levels of triggering 
were used (see rcf Q for a detailed description). The 
spectrometer included two fan-shaped arrays of scintil- 
lator hodoscopes which provided fast triggering of the 
first level by requiring a minimum transverse momen- 
tum. This trigger passed a typical event rate of 100 KHz 
with a transverse momentum cut-off resolution of about 
7%. The second level trigger selected high transverse mo- 
mentum particles with a 4% momentum resolution and 
accepted a rate of 10 KHz. The third level trigger allowed 
a wide range of choices for the angles of the two fast par- 
ticles. The accepted rate ranged from 10 to 40 /sec. 

The coordinate system was chosen with the z coor- 
dinate in the beam direction and the y direction nor- 
mal to the scattering plane (x, z). The latter is defined 
by the incident beam and either one of the emerging 
protons. The data were analyzed in terms of the mo- 
menta in the y direction {Pf y ) and the light cone variable 
a = (Ef — Pf z )/m, where Ef is the total target nucleon 
energy, Pf z is the z component of the target proton mo- 
mentum and m is the mass of the target proton. The 
variable a is a natural choice for high energy reactions 
and is also ideally adapted to our experimental analysis. 
We determined a with a precision of a ~ 3%. Setting 
Ef ~ m, we can write a ~ 1 — Pf z /m, and with the 
additional good approximation in our kinematical region 
Eq = Pq, we can write: 

s ~ to 2 + to 2 + IraP^a (1) 

where Eq and Pq are the energy and momentum of the 
incident proton and ifi is the off-mass-shell mass of the 
target nucleon. Setting m ~ m: 

s ~ 2m 2 + 2mP a (2) 

This makes s proportional to a. The Pf y had a res- 
olution of a = 40 MeV/c and the resolution in Pf x was 



(7 = 170 MeV/c. Because of its better resolution, Pf y 
was used to represent a transverse component. Both the 
a and Pf y resolutions were determined by the elastic pp 
scattering events from the CH 2 target. 
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FIG. 1. Emiss distribution for 6 GeV/c, 0.1 < \P fy \ < 0.2 
GeV/c, and |«-0.87| < 0.05. The Gaussian (dotted line) rep- 
resents the QE events. The shape of the background (dashed 
curve) was deduced from extra track events. The fall-off above 
1.5 GeV is an artifact of the trigger. The solid curve is the 
result of the fit (see text). 

The quasi-elastic scattering nature was ensured by se- 
lecting events from the C target with just two tracks in 
the detector and by software cuts on the quality of track 
reconstruction. An upper limit on the excitation energy 
of the residual nucleus ( E m i SS ) was imposed in order 
to eliminate events where additional particles could be 
produced. Given our resolution for E m i SS , we applied a 
cut of |-E m i SS | < 0.5 GeV. Since this cut is above m^, 
some inelastic background from soft neutral particles, 
such as those coming from pA — > ppir (A — 1) events, 
could penetrate the cuts and had to be subtracted. The 
shape of this background was determined from a fit to 
the E m i ss distribution of events with extra tracks in the 
spectrometer. Subsequently, we used that shape to fit 
the E m i SS spectrum for 50 MeV/c wide Pf y bins to the 
background (BG) and a Gaussian centered at E m i SS ~ 
which represents the QE events. The width of the Gaus- 
sian was determined from a fit to the peak at E m i SS ~ 
for < Pf y < 50 MeV/c, where the peak is very promi- 
nent. The reported a and Pf y distributions are the re- 
sults of these fits. An example is shown in Figure [l]. 
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FIG. 2. a-b: |P/ S | distributions for different a slices. The 
6 GeV/c (solid line) and the 7.5 GeV/c (dashed line) distri- 
butions were normalized at Pf y — 0. c: \Pf y \ distributions 
for a = 0.79 (solid line) and for a = 0.87 (dashed line). The 
data in each a range are summed for 6 and 7.5 GeV/c (see 
text). The two distributions were normalized at Pf y = 0. 

In Figure || we present \Pf y \ distributions. Figs. ||(a)- 
§(b) show the Pf y distributions for two regions of a, at 6 
and 7.5 GeV/c. The curves were normalized at Pt y — 
and we see that the shapes are independent of the inci- 
dent energy within the experimental errors. This is what 
one expects from the impulse approximation (IA). In the 
IA the QE cross section is factorized into contributions 
from the probability to find a nucleon with Fermi mo- 
mentum Pf in the nucleus n(Pf z , Pft), the free pp elastic 
cross section and the nuclear transparency T(s,t): 



d 3 a 
dP f 



(s,t) q 



da 



(s,t)f ree ■ n(P fz ,Pf t ) ■ T(s,t)dt 



(3) 



For scattering at 9 cm ~ 90° the free cross section is a 
function of s only. For a narrow region of a corresponding 
to a narrow region of s, we can then write: 



d 2 a 
dP}t 



(«)<?. e. = n(Pft) ■ Factor(a) 



(4) 



This means that for fixed a ,the \Pf y \ distributions at 
the two incident energies will scale by some function of s 



as can be seen clearly in Figure || ■ Note that for small 
a (Fig. |(b)) the P fz is about 200 MeV/c and, since 
the scaling holds also for large Pf y , the scaling has been 
checked up to a fairly large Fermi momentum (about 0.5 
GeV/c). 

In Fig. ^|(c) we compare the transverse Fermi dis- 
tributions for two different a ranges. Since the shapes 
were shown to be independent of the incident energies 
we summed the measured distributions for the two in- 
cident energies. We observe that the transverse Fermi 
distribution gets wider for larger longitudinal Fermi mo- 
mentum distributions (smaller a). This is in contradic- 
tion to what one would expect from a simple Fermi gas 
model. With a given distribution \Pf\, a large Pf z would 
give a narrower Pf y distribution. 

As we mentioned in the introduction, the large momen- 
tum transfer QE (p, 2p) reaction prefers small s (small 
a) due to the strong s dependence of the elementary pp 
elastic cross section. The flux factor for the quasi-elastic 
pp cross section is different from the free pp one, due to 
the motion of the target proton in the nucleus. Since 
the differential cross section is not invariant, there are 
also Jacobians that have to be included when changing 
reference frames. Taking these effects into account, we 
obtained the nuclear momentum distribution, n(Pf z ), by 
multiplying the measured a distribution, -4^, by the fac- 



tor of Equation |5[ 
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where ^ pp is the measured free pp cross section [ pT| , 
Pi is the momentum of one of the outgoing protons, Ei 
and Pi z are the energy and z momentum component of 
the other one, and Pit and P2t are the transverse compo- 
nents of the two outgoing particles. A derivation of the 
factor is given elsewhere || . 

After this correction the two measured a distributions 
at 6 and 7.5 GeV/c were consistent with each other in 
shape, so we joined the two sets of measurements. The 
result is shown in Figure ^ which, up to distortion from 
initial and final state interactions, represents the Pf z dis- 
tribution of the target nucleon. 

The transverse distribution (Pf y ), shown in Figure]^, 
was obtained by adding the measured distributions of 
the measured a regions, weighted according to the inte- 
grated number of events in the measured a distribution. 
We see that the longitudinal and the transverse distribu- 
tions have the same shape in spite of the very different 
procedures used to get them. 

Also in Figure ^, the data are compared to a prediction 
of a simple independent particle Fermi motion distribu- 
tion. An harmonic oscillator model (HO) was used with 
parameters fitted to 12 C from electron scattering with 
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the kinematical constraints of our measurement p2[ . As 
can be seen clearly, the independent particle model fails 
substantially to describe the large momentum tails of 
the distribution. In order to quantify this statement, 
we present the ratio between the number of events with 
154 < Pf z < 280 MeV/c and the number of events with 
< Pf z < 154 MeV/c (we measured events with Pf z > 
and we assumed that the distribution is symmetric about 
P fz = 0). The measured ratio is (29 ± 2)%. The HO 
prediction for the same ratio is only 11.7%. 
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FIG. 3. Pfz is the longitudinal momentum distribution, ob- 
tained from the a distributions measured at 6 and 7.5 GeV/c 
and corrected for the s dependence induced by the elemen- 
tary free cross section. \Pj y | is the transverse distribution, ob- 
tained from the measured transverse distributions, at several 
a regions (see text) . HO is a harmonic oscillator independent 
particle model calculation. The Pf z and HO distributions are 
normalized to 1000 at the first bin. 

In conclusion, the QE events were identified and we 
showed that it is possible to separate the nuclear prop- 
erties (Fermi momentum distribution) from the reaction 
mechanism up to a total Fermi momentum of about 0.5 
GeV/c. Based on this result, we deduced transverse and 
longitudinal Fermi distributions up to large momenta. 
The tails of these distributions are consistent with elec- 
tron scattering data Jl|| and are larger than predicted 
by independent particle models and are characteristic for 
short range correlation. We also found that the trans- 
verse Fermi distribution increases as Pf z increases. 
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